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A comparison is made between results obtained using smooth initial 
conditions and event-by-event initial conditions in the hydrodynamical de- 
scription of relativistic nuclear collisions. Some new results on directed flow 
are also included. 
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1. Objective 

Hydrodynamics has been rather successful at describing data obtained in 
relativistic nuclear collisions at RHIC. Usually, smooth initial conditions are 
assumed (see e.g. fig.l in pQ and fig. 3 in [2J. On the other side, microscopic 
codes such as NeXus predict initial conditions event-by-event, which are 
quite irregular as shown in fig.l. 

The question we address here is whether such structures (hot spots or 
more precisely hot tubes) can have a sizable effect on variables. 

To solve the hydro equations with very irregular initial conditions, we 
use the SPheRIO code. This code is based on the method of Smoothed 
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Au+Au, evento aleatorio Cu+Cu, evento aleatorio 




Fig. 1. T] = slice for initial energy density of a RHIC collision in the 6-15 % 
centrality window. 

Particle Hydrodynamics, originally developed in astrophysics and adapted 
to relativistic heavy ion collisions in [3] . The version of NeXSPheRIO used 
here has initial conditions provided by NeXus [1] and normalized by an 
77-dependent factor to reproduce dN^/drj in each centrality window [5]. 
The equation of state has a critical point [6j. Tj out is fixed (mostly) by 
dN c h/p t dpt and depends on the centrality window (i.e. number of partici- 
pants). Centrality windows are defined using participant number and not 
impact parameter [7J. An ideal fluid is assumed, a code with Smoothed 
Particle Hydrodynamics and dissipation is under development [8]. 

2. Comparison between fluctuating and average IC 

In the following, we present a summary of results obtained using smooth 
initial conditions and running once the SPheRIO hydro code (standard ap- 
proach) or using a set of NeXuS initial conditions, running for each initial 
conditions the SPheRIO hydro code and computing averages over the set 
for various observables (event-by-event hydrodynamics). 

2.1. pt distribution 

As can be seen in figure 2 (left), the high pt part is lifted. This is 
expected since hot tubes must expand more violently, producing more high 
Pt particles [9j [10] . 

2.2. elliptic flow 

V2(pt) is flatter as seen in figure 2 (centre). This is also expected as the 
isotropic expansion of hot tubes produces more high pt particles and lowers 
vzipt) 13 HD]- In addition, ^(f?) has no shoulder [11] as seen in figure 2 
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Fig. 2. Left: charged particle pt distribution. Solid line: e-by-e initial conditions. 
Dashed: smooth initial conditions. Data: [T^]. Center: p t dependence of (t^). 
Data: [T3]. Right: n dependence of (t^)- Data: [T5] , 

(right). The effects (isentropic expansion) of the hot tubes are more visible 
in regions of lower matter density present at larger n's I10j . 

2.3. Other comparisons 

In [14] . we argued that the hot tubes should manifest themselves giving 
smaller HBT radii. However, the situation might be more complicated. 

Another observable where hot tubes might manifest themselves is the 
ridge, a structure observed in the 2 particle correlations, plotted as function 
of pseudorapitity difference Ary and azimutal angle difference A(p between 
a high pt trigger hadron and its associated hadrons (see e.g. [IS])- The 
structure is An independent. In NeXSPheRIO, the hot tubes can lead to 
such a ridge for the e-by-e initial conditions and not the smooth ones |16j . 

Finally, the fluctuations in the e-b-e initial conditions also manifest 
themselves in fluctuations of V2 (as well as v±). The predicted values for 
V2 at 130 A GeV |17j and estimates at 200 A GeV [5] are in agreement 
with data |18[ I19j . Improvements to remove the non-flow effects have been 
reported by STAR and PHOBOS, see e.g. [20]. 

3. New results on directed flow 

In this section, we present some new preliminary results obtained with 
NeXSPheRIO on directed flow. 

3.1. What is directed flow and what is expected 

If a nucleus-nucleus collision is a number of independent nucleon-nucleon 
collisions, the momentum distribution is isotropic. If instead, it leads to 
thermalized matter in the overlap region, the momentum distribution is 
stretched along the impact parameter direction, V2 is a measure of this 
stretching (so teaches about IC, thermalization, etc). There is also the 
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possibility that the momentum distribution be shifted/deformed towards 
one of the sides in the x-y plane, V\ is a measure of this shift. 

At some energy, a "wiggle" in v\{rj) is predicted. In some microscopical 
models such as RQMD and UrQMD, this could be the case for nucleons at 
RHIC energy [2TJ[22]. In hydro models, this could be the case for the fluid, 
if a QGP phase occurs [231 121 US E£5] • 

At SPS energy (40 A GeV and 158 A GeV), it was shown [27] that pions 
and protons behave oppositely. Pion directed flow as function of rapidity 
has no wiggle and crosses y=0 with a negative slope while nucleon directed 
flow has no wiggle and crosses y=0 with a positive slope (except perhaps 
at the higher energy, in the more peripheral bin, where there is a hint of 
wiggle). 



3. 2. RHIC results on directed flow 

At RHIC, directed flow for charged particles is rather similar to what was 
obtained at SPS for pions: it crosses r/ = with a negative slope [28 ], [29] , [30] . 
This is understandable since charged particles are mostly pions, the fluid 
directed flow must be dominated by pions. The turnover in v\ (77) occurs for 
different values of r\ in PHOBOS and in STAR (see below). 

Results for identified particles are becoming available [31j . 

In addition, comparison of results for directed flow in Cu+Cu and Au+Au 
collisions show no system-size dependence [30] . 



3. 3. NeXSPheRIO results on directed flow 

NeXSPheRIO results are in qualitative agreement with PHOBOS for 
all 77's and quantitative agreement for | 77 |< 3 (figure 3 left). They are 
in qualitative agreement with STAR for | 77 |< 3 but turnover occurs for 
smaller 77 than for STAR (figure 3 right). 



v,> for AuAu at s =200 GeV, 0-40 % 



1 PHOBOS, Back etui 
■ NeXSPheRIO 3954 events 











Fig. 3. Comparison of charged particle (vi) for NeXSPheRIO with (left) PHOBOS 
[28] and (right) STAR [30]. 
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v i (rj) from NeXSpheRIO for various centrality windows for Au+Au and 
Cu+Cu at 200 A GeV is shown in figure 4. Little dependence on A is seen 
in the windows 6-f 5% to 45-55%. Statistics must be improved. 




Fig. 4. Comparison of (vi) obtained in Cu+Cu and Au+Au, from NeXSPheRIO. 

Figure 5 (left) illustrates particle dependence. In NeXSPheRIO, protons 
have a big wiggle, pions have a plateau (left). A similar result was obtained 
using UrQMD [22]. In figure 5 (right), it is seen that vi(r]) has a plateau for 
fluctuating initial conditions and a somewhat stronger negative inclination 
for smooth initial conditions. 




Fig. 5. Left: (vl) for pions and protons. Right: (vi) for e-by-e and smooth initial 
conditions. 



4. Summary 

A short review of possible effects of fluctuating initial conditions, rather 
than smooth ones, was presented. In addition to providing a reasonable 
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description of various observables, as is possible with smooth initial condi- 
tions, some new effects were listed, most notably the ridge effect and the V2 
fluctuations, which do not appear when using the smooth initial conditions. 
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